An efficient polynomial time approximation scheme for load 
balancing on uniformly related machines 

Leah Epstein* Asaf Levin ^ 

Abstract 

We consider basic problems of non-preemptive scheduling on uniformly related machines. For 
a given schedule, defined by a partition of the jobs into m subsets corresponding to the m machines, 
PsJ I Ci denotes the completion time of machine i. Our goal is to find a schedule which minimizes or 

^— s \ maximizes ^ Cf for a fixed value of p such that < p < oo. For p > 1 the minimization problem 

. is equivalent to the well-known problem of minimizing the Ip norm of the vector of the completion 

' times of the machines, and for < p < 1 the maximization problem is of interest. Our main result 

^ . is an efficient polynomial time approximation scheme (EPTAS) for each one of these problems. 

Our schemes use a non-standard application of the so-called shifting technique. We focus on the 
OO . work (total size of jobs) assigned to each machine and introduce intervals of forbidden work. These 

intervals are defined so that the resulting effect on the goal function is sufficiently small. This allows 
the partition of the problem into sub-problems (with subsets of machines and jobs) whose solutions 
C/^ ■ are combined into the final solution using dynamic programming. Our results are the first EPTAS's 

. for this natural class of load balancing problems. 

O 

1 Introduction 

J> I We consider non-preemptive scheduling problems on m uniformly related machines. In such problems, 

^ ■ we are given a set of jobs {1,2,..., n}, where each job j has a positive size pj. The jobs need to be 

^ , partitioned into m subsets Si, ... , Sm, with Si being the subset of jobs assigned to machine i. We let 

t:J" \ Si denote the speed of machine i, and the processing of a job j takes ^ time units if j is assigned to 

^ ■ machine i. For such a solution (also known as a schedule), we let Ci = be the completion 

. tinie of machine i. The work (also called the weight) of machine z is VFj = Yl,j^s Pi ~ ^« ' 

that is, the total size of the jobs which are assigned to i. The makespan of the schedule is maxj Ci, 
and the optimization problem of finding a schedule which minimizes the makespan is well-studied (see 



^ . e.g. [ 20l[T9lE3l|24l|26l). The problem of finding a schedule which maximizes mini Ci is the well-known 

\ Santa Claus problem on uniformly related machines (see e.g. |[T8l[30l l2ll5l[T5l[8l). Both these problems 

^ ■ are concerned with the optimization of the extremum values of the set {Ci, . . . , Cm}- 

Motivated by minimizing average latency in storage allocation applications (rather than worst-case 
latency), researchers have suggested to study the optimization goal of minimizing the £2 norm (and the 
goal of minimizing the ip norm for p > 1) of the vector of completion times of the machines (see 
e.g. |[T2l[mi28l l4l[3l). It was stated more recently by Bansal and Pruhs [Tl that: ''The standard way to 
compromise between optimizing for the average and optimizing for the worst case is to optimize the ip 
norm, generally for something like p = 2 or p = 3." An additional perspective of using the £p norm as 
an objective function has arisen recently in algorithmic game theory [9]. Note that the minimization of 
the £p norm is equivalent to minimizing the sum of the p-th powers of the completion times of machines. 
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Thus, we consider objective functions in which the entire vector C = (Ci, . . . Cm) affects the value of 
the objective function. Our class of objective functions includes the minimization of the sum of the p-th 
powers of the completion times of machines which is equivalent to the minimization of the £p norm of 
C. More precisely, given a fixed real (finite) number p such that < p < oo, we consider the problem of 
minimizing Yl^i ^^'^ problem of maximizing X^^^ Cf. The minimization problem for p < 1 
is trivially solved by placing all the jobs on one of the fastest machines. Therefore, we consider the 
minimization problem only for values of p such that p > 1. Similarly, the maximization problem is 
trivially solved for p > 1 by placing all the jobs on one of the slowest machines. Hence, we consider 
the maximization problem only for values of p such that p < 1. 

An 7?--approximation algorithm for a minimization problem is a polynomial time algorithm which 
always finds a feasible solution of cost at most TZ times the cost of an optimal solution. An TZ- 
approximation algorithm for a maximization problem is a polynomial time algorithm which always finds 
a feasible solution of value at least ^ times the value of an optimal solution (we use the convention of 
approximation ratios greater than 1 for maximization problems). The infimum value of TZ for which an 
algorithm is an 7^-approximation is called the approximation ratio or the performance guarantee of the 
algorithm. A polynomial time approximation scheme (PTAS) is a family of approximation algorithms 
such that the family has a (1 + e) -approximation algorithm for any s > 0. An efficient polynomial 
time approximation scheme (EPTAS) is a PTAS whose time complexity is of the form /(-) • poly{n) 
where / is some (not necessarily polynomial) function and poly{n) is a polynomial of the length of the 
(binary) encoding of the input. Motivated by this definition of polynomial time complexity, we say that 
an algorithm (for some problem) has polynomial time complexity if its time complexity is of the form 
f{-) ■ poly{n). Note that whereas a PTAS may have time complexity of the form n^^~\ where g is for 
example linear or even exponential, this cannot be the case for an EPTAS. The notion of an EPTAS is 
modern and find its roots in the FPT (fixed parameter tractable) literature (see [ 10llT3llT7ll29l ). 

Our main result is a class of EPTAS's for minimizing Yl^i ^'^^ fixed value ofp> 1, and for 
the problem of maximizing YliLi Cf for any fixed positive value of p < 1. Note that these problems 
are known to be strongly NP-hard even for identical machines (via the standard reduction from the 3- 
PARTITION problem) and therefore our results are the best possible. Our results are the first EPTAS 's 
for these important load balancing problems on uniformly related machines. 

The running time of an EPTAS (and of a PTAS) is expected to be polynomial in the number of 
jobs as well as in the number of machines. For a fixed (constant) number of machines, load balancing 
problems typically have a fully polynomial time approximation schemes (EPTAS's, which are EPTAS's 
where / is polynomial) ll25ll6l[T4l[T6ll. 

We next review the previous PTAS and EPTAS results for an arbitrary (non-constant) number of 
uniformly related machines and the special case of identical machines (where all machines have unit 
speed). It was shown by Hochbaum and Shmoys that the makespan minimization problem has a PTAS 
for identical machines ||23T| and for uniformly related machines [24 1. It was noted in ||2D that the 
PTAS of Ii23l for identical machines can be converted into an EPTAS by using integer program in fixed 
dimension instead of dynamic programming. Recently, Jansen (26] was able to solve the long-standing 
open problem of establishing an EPTAS for the makespan minimization problem on uniformly related 
machines. The Santa Claus problem is also known to have a PTAS and an EPTAS for identical machines 
Il30l l2l. For uniformly related machines a PTAS is known ll5l[T5l. 

The problems studied here are known to have an EPTAS on identical machines fP, "21, and a PTAS 
on uniformly related machines 1 15 1. The existence of an EPTAS for these problems on uniformly related 
machines was stated as an open problem by lITSl . This open problem is resolved in our work. 
Outline. Our EPTAS's have the following structure. First, we sort the machines in a non-decreasing 
order of their weight in an optimal solution (according to either non-increasing or non-decreasing speed). 
We note that some machines may get a zero weight; we guess their number and remove those machines 
from the instance. We round the processing times of the jobs and the speeds of the machines, so that 
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the number of possible values is reduced sufficiently, and so that all job sizes are integer multiples of a 
small value. 

Next, we observe that we can extend the EPTAS for identical machines to the case where we are 
guaranteed that in an optimal solution the ratio between the maximum work of any machine and the 
minimum work of any machine is bounded. We show that in this case the speed ratio is bounded as well. 
We extend this EPTAS further to allow some total size of jobs to remain unscheduled. This will be our 
building block in the design of the EPTAS for the general case. 

To reduce the general case into a series of sub-problems of the former type, we create gaps between 
the set of allowed weights of machines. For that, we apply the so-called shifting technique |[22l in an 
original way. Afterwards, we apply dynamic programming to determine the series of sub-problems, that 
is, the intervals of machines whose weights come from each interval of allowed weights. The EPTAS 
for the special case is used as a black box in this dynamic programming, where unscheduled jobs of one 
sub-problem are scheduled later by another sub-problem. Omitted proofs can be found in the Appendix. 

2 Preliminaries 

In this paper we consider the sum of the p-ih powers of a vector rather than the (|)-th power of this 
value. Note that since p is a fixed constant, our results apply also for this last alternative measure (which 
is the ip norm for the case p > 1). Throughout the paper, for a solution A we denote by A both the 
solution and the value of the objective function for this solution. 

When we consider the maximization problem, we sometimes allow the algorithm to avoid assigning 
some of the jobs. It is clear that adding these jobs arbitrarily to the schedule can only improve the 
solution. Hence, if we can bound the total value of the solution which assigns a subset of the jobs, after 
adding the unscheduled jobs (to create a complete solution), we get (at least) the same performance 
guarantee. 

Let e be a small constant such that < e < ^ and ^ is an integer. Epstein and Sgall [.15.1 observed 
the following claim. 

Claim 1 Let ii and 12 be a pair of machines such that si^ < Si^, that is, 12 is faster than ii. Consider 
the minimization problem when p > 1, then any optimal solution satisfies Wi-^ < Wi^. Consider the 
maximization problem when p <\, then any optimal solution satisfies Wi-^ > Wi^. 

Motivated by the above claim we will sort the machines according to their weights. That is, when 
we consider the minimization problem we will assume that si < S2 < ■ ■ ■ < Sm, whereas when we 
consider the maximization problem we will assume that si > S2 > ■ ■ ■ > Sm. In this way, machines of 
lower indices should get smaller weight than machines with higher indices (or equal weight). We next 
consider a pair of machines ii and 12 such that is significantly smaller than Si^. We know that in the 
minimization problem Wi-^ < Wi^ and in the maximization problem Wi^ > Wi^ . Our next goal is to 
strengthen these bounds. Let 5 be such that Q < 5 < e. 

Lemma 2 Consider the minimization problem (p > 1), and a pair of machines i < i'. There exists a 
function a{6) = | such that if si < a{5) ■ Si' then in any optimal solution Wi < 6Wi'. 

Lemma 3 Consider the maximization problem (p < 1), and a pair of machines i < i'. There exists a 
function a{5) = ((1 + 5Y — 1)^^^ such that ifsi' < a{6) ■ Sj then in any optimal solution Wi < 6Wii. 

Note that a{6) < 6 < e. This is clear for the minimization problem, and for the maximization 
problem it holds because (1 + 6y < 1^ + 6^ = 6^ + 1 where the inequality holds by the concavity of 
xP for p < 1, and the claim holds by the monotonicity of x^^p. We summarize the last two lemmas by 
the following straightforward corollary, which we will use. 
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Corollary 4 Consider a pair of machines i and i' such that i<i'.IfWi> ^^Wiifor some function 7 
(such that 7(e) > ^ ), then there is a function P such that the ratio between the speeds max{sj, Sj/} and 
min{si, Sj/} is bounded by /3(e). 

First rounding step. In what follows we assume without loss of generality that the speeds are integer 
powers of 1 +e. This assumption is justified by the observation that increasing the speed of each machine 
to the next value of the form (1 + e)-' (for integer j) may decrease the completion time of this machine 
by a multiplicative factor of at most 1 + e. Thus approximating the optimization problem with respect 
to the new speeds within a factor of 1 + e gives a (1 + e)^"*"^ approximation to the original instance of 
the problem. Thus by scaling e accordingly, the assumption is justified. Moreover, we assume also that 
the sizes of all jobs are integer powers of 1 + e. This assumption is justified by the observation that 
increasing the size of each job to the next value of the form (1 + ey (for an integer j) may increase the 
completion time of each machine by a multiplicative factor of at most 1 + e and may not decrease it. 
Thus the following properties can be assumed. 

Assumption 5 The speed of each machine as well as the size of each job is an integer power ofl + e. 

Second rounding step. Letpmax = maxj=i 2,...,nPj. We will let OPT denote a fixed optimal solu- 
tion for the resulting instance after the second and final rounding step which we now define. For the 
minimization problem, we apply the following rounding (down) of the (rounded) processing times. If 
Pj < 22sax ^gjj round pj down to be zero and we remove all such jobs from the instance. Otherwise, 

2 

we round pj down to the next integer multiple of fj, = ^ Given a solution for the input after this 
second rounding step, we create a solution for the original instance, by assigning all the removed jobs to 
the machine in which a job of size Pmax is assigned (breaking ties arbitrarily). 

Lemma 6 The cost of every solution to the instance of the minimization problem after the second round- 
ing step is no larger than its cost before this rounding step. Moreover, given a solution to the final in- 
stance, the cost of the resulting solution for the input after the first rounding step is at most (1 + e)^^ 
times its cost for the new instance. 

Proof. Since we only round down processing times of jobs, the cost of the new solution cannot increase. 
Therefore, the first claim holds. As for the second claim we first consider the effect of returning the 
removed jobs to the solution. The work of the machine which receives these jobs increases by at most 
epmax and since it was previously at least Pmax, it increases by a multiplicative factor of at most l + e. 

Next, consider the effect of reverting the processing time of the jobs whose processing times were 
larger than SSsss ^Yiqit values in the instance after the first rounding step. Note that such a processing 

2 

time of a job is increased by an additive factor of at most ^ ^^'^'^ and hence by a multiplicative factor of 
at most l + e. Therefore, the work of each machine is increased by a multipHcative factor of at most 
l + e, and the claim follows. ■ 

For the maximization problem, the second rounding step is defined as follows. In this case we let 

1 

= ^^'^"f"' , and we round the processing time of every job down to an integer multiple of p. 

n-mP 

Lemma 7 The value of every solution to the new instance of the maximization problem after the second 
rounding step is not larger than its value before this rounding step and thus reverting jobs to their sizes 
after the first rounding step can only improve the performance. Moreover, consider an optimal solution 
SOL to the instance after the first rounding step. Denote by SOL„e„,., SOLo/^^ its objective fiinction values 
in the new instance after the second rounding step, and before the second rounding step, respectively. 
Then, OPT > SOL„e«, > (1 - e)SOLo,d. 
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Note that after the second rounding step the size of any job is an integer multiple of /i. 

Lemma 8 Given an interval [L, U], the number of distinct sizes of jobs is at most logi_|-£ + 2. 

Proof. At the end of the first rounding step the number of distinct sizes of jobs in this interval is at most 
log^^g + and this number may increase by at most 1 due to the second rounding step (in case where 
we round down a value which is shghtly larger than U). ■ 

3 Approximating the problem with a bounded weight ratio 

In this section we consider the following variant of our (maximization or minimization) problem which 
is called BOUNDED RATIO (BR). The input to this problem consists of the following parts: 

1. A set of £ > 1 consecutive machines with speeds Sj, Sj+i, . . . , for which we define 

Smin = min{si, Sj+i, . . . , Sj+£_i} and Smax = max{si, Sj+i, . . . , Sj+^-i} , 

and assume < /3(e). Recall that for the minimization problem we assume that the speeds are non- 
decreasing and thus Smin = Si and Smax = Si^i-i, and for the maximization problem we assume that 
the speeds are non-increasing and thus Smin = and Smax = 

2. A pair of values fi < Wi < yVi+e-i bounding the weights of machine i and machine i + £ — 1, 
respectively, such that ^yj~^ < 7(e)- Wj, Wj+^^i are (not necessarily positive) integer powers of 
1 + e. 

3. A value A which is an integer multiple of /i. For the minimization problem A is an upper bound on 
the total size of jobs that the algorithm does not necessarily need to assign to any of these machines, 
whereas for the maximization problem the value of A is a lower bound of the total size of jobs that the 
algorithm should not assign to any of these machines. 

4. A set L of large jobs 1, 2, . . . , n each of size at least eWi. 

5. A value B of the total size of existing small jobs . The small jobs can be assigned fractionally. The 
value B is an integer multiple of fi. 

The goal is to schedule the large jobs and the small jobs on the £ machines such that the weight of 
each machine is at least Wi and at most Wj+^-i. We allow an arbitrary subset of the jobs (out of the 
large jobs and small jobs) to remain unscheduled as long as in the minimization problem its total size is 
at most A, and in the maximization problem the total size of the unscheduled jobs must be at least A. 
We assume that such an assignment of the jobs is feasible or else the algorithm returns FALSE. The goal 
is to minimize or maximize the value Cj of the schedule, and the input to BR consists of i, £, 

Wi, Wi+e~i, A, L and B. We will allow the algorithm (not the optimal solution) to use machines with 
weights in the interval i^^-, Wj+^-i • (1 + 3e)]. 

Remark 9 By Corollary^ the requirement that < R(e) follows from ^'+f~^ < 7(e). 

^min ''^i 

Next, the number of different sizes of large jobs is at most logi+j + 2 < log^^^ + 2 and 

thus this number is a function of e. Let H be the set of different sizes of large jobs, and for each h ^ H 
we let Uh be the number of jobs of size h. We define a class of machines to be machines with the same 
(rounded) speed. Since the speeds are integer powers of 1 + e and the ratio between speeds satisfies 
^ /3(e) we conclude that the number of non-empty machine classes, denoted as r(e), is at most 

^min 

log^^j /3(e) + 1 that is a function of e. We denote the non-empty machine classes in our problem by 
Ml, . . . , -Mr(g). For each machine class M^, whose machines have a common speed of au, we denote 
by vic^k) = l-^fcl the number of machines in Mj.. 

We define a configuration K of a machine as a vector with the following components. The first 
\H\ components of K define the number of large jobs of each size which we schedule on a machine 
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with configuration K. For each h £ H, n{h, K) denotes the number of jobs of size h which are 
scheduled on a machine with this configuration. Each K) is a non-negative integer which is at most 
^ly^T^ < , that is a function of e. The next component of K is an integer power of 1 + e in the range 
[Wj, Wt+£_i] denoted as w{K). For the minimization problem, let w{K) = [ "'^^^ J • /x. w{K) is the 

maximum total size of jobs in configuration K. For the maximization problem, let w{K) = [ ^^^-* ] • ^u. 
w{K) is the minimum total size of jobs in configuration K. The number of options for this component 
is at most logi+g ^y^^~^ + 1 < ^^^i+s + 1' that is a function of e. The last component is the 
machine speed, and we denote this component by s{K). There are r(e) options for this last component. 
We conclude that the number of different configurations is a function of e and we can enumerate all of 
them in a constant time. We denote by K, the set of all configurations. A configuration of a machine 
defines the number of large jobs of each size which are scheduled on such a machine, as well as the total 
size of small jobs which are scheduled fractionally on such a machine (which is the difference between 
w{K) and the total size of the large jobs). Since the size of a job is an integer multiple of ji, we require 
that the total size of small jobs which are scheduled on such a machine is an integer multiple of ji as 
well. 

We define an integer program of fixed dimension to solve BR. The decision variables are for each 
configuration A' G /C, a variable xk counting the number of machines which are scheduled according 
to configuration K. We let yh be the number of large jobs of size h £ H which remain unscheduled in 
our solution. 

The following integer program is used for solving our minimization problem. 

s.t. J2KeK:s{K)=a, = ^i^k) V/c = 1, 2, . . . r(e) (1) 

YjKeK "-(^i K) ■ XK + Vh = n-h e H (2) 

E/je// h-yh- J2k€K ■ {w{K) - Y.h&H ^ ' ^i^^ K)) <A-B (3) 

Y.heHh-yh<A (4) 

XK,yh>^ yK eicyhe H. (5) 

The family of constraints ([T|l enforce that we use only v{(Jk) machines with speed ak- The family of 
constraints ^ enforce that exactly yh jobs of size h are unscheduled by our solution. Constraints Q 
and dll) enforce the condition on the total size of jobs which are unscheduled. To see this last claim first 
note that a machine which is scheduled with configuration K leaves a gap of size w{K) — J2h&H ^ ' 
n{h, K) for a possible scheduling of small jobs (since all jobs have sizes which are integer multiple of 
^). There are two cases. In the first case the total size of the gaps (of all the machines) is sufficient 
for scheduling all the small jobs that is YIk^k. 

case, 

we assume without loss of generality that all the small jobs are scheduled. Thus in this case we only 
need to make sure that constraint ^ holds. This last constraint holds because the total size of the 
large jobs which are unscheduled is exactly J2heH ^ ' Vh^ ^^'^ thus such a feasible solution satisfies 
both constraints ^ and (IDl. In the other case there are small jobs which are not scheduled by the 
solution. Since we allow fractional scheduling of small jobs, their total size is exactly B — YIkgJC ' 
(^w{K) — J2heH ^ ■ "^(^1 ^)) which is positive. Thus the total size of the large jobs which may be 
unscheduled is at most A — — J2k£K ' (^(-^) ~ "l^h^H ^ ' ^(^' ^))) ^^'^ '■h^^ constraint ^ 
holds, and by the assumption of this case, we conclude that constraint ^ holds as well. 

Let (x* , y * ) be an optimal solution for the above integer program, and let X* = Yl ( "^(k) ) ' 
be its objective function value. 

Claim 10 Denote by OPT^r the optimal solution for this bounded ratio minimization problem as well as 
the value of its objective function. We have X* < (1 + eY ■ OPTfo^. 
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The following integer program is used for solving our maximization problem. 

max 

s.t. T.KeK:s(K)=ak^K = v{crk) VA; = 1, 2, . . . r(e)(6) 

YjKeK "'(^' K) ■ XK + yh = rih yh£ H (7) 

EheH h-yh-Zxe^XK- {w{K) - ^^^^ h ■ n{h, K)) > A - B (8) 

T.KeK ■ {w{K) - ZheH h ■ n{h, K)) < B (9) 

XK, yh>0 WK elC,yhe H.10) 

The family of constraints ^ enforce that we use only v{ak) machines with speed cjf^. The family 
of constraints ([7]) enforce that exactly jobs of size h are unscheduled by our solution. Constraints 
dH) and ^ enforce the condition on the total size of jobs which are unscheduled. To see this last 
claim first we observe that a machine which is scheduled with configuration K leaves a gap of size 
w{K) — J2heH ^ ' "^(^' ^) which must be covered by small jobs. Therefore, the total size of small jobs 
should be sufficient to fill all these gaps. This enforces constraint Q. The remaining small jobs together 
with the unscheduled large jobs need to be of total size of at least A. Thus constraint ([8]l holds as well 

Let (x* , y*) be an optimal solution for the above integer program, and let X* = ( "^(x) ) 
be its objective function value. 
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Claim 11 Denote by OPT;,,, the optimal solution for this bounded ratio optimization problem as well as 
the value of its objective function. For the maximization problem X* > 



OPTb, 



For each of our problems, we first solve the corresponding integer program. We show that given the 
solution {x* ,y*) to the integer program, we can construct a feasible solution for BR whose objective 
function value is at least as good as X* . Then large jobs are assigned to the machines according to the 
configurations of the machines. Small jobs are assigned using (fractional) next-fit to the remaining gaps. 
In fractional next-fit we assign the jobs one by one until the current gap does not have a sufficient room 
for the next job, in which case we assign a fraction of the job to the current gap, so as to fill exactly the 
gap, and the remaining of the current job is assigned to the next gap (of the next machine). The work 
of a configuration K is assumed to be w{K). This procedure fills exactly all the gaps until one of the 
following two cases occurs. Either, there are no additional machines and there are still small jobs to be 
assigned, or there are no additional small jobs, but remaining gaps. By constraint I© in the maximization 
problem only the second case may occur. In the minimization problem the first case does not cause any 
problem as it gives a feasible solution to problem BR since the remaining small jobs are counted towards 
A in Constraint ([3]). In the second case for maximization problem, the unscheduled small jobs (or parts 
of these) are counted towai^ds the total size of unscheduled jobs as well. This assignment algorithm 
is clearly a polynomial time algorithm. By the above claims it suffices to show that we can solve the 
integer programs in polynomial time. 

Claim 12 The integer programs can be solved in strongly polynomial time. 

Proposition 13 Problem BR has an EPTAS. 

We next consider a variation of problem BR in which the small jobs which are scheduled on one 
of the £ machines need to be scheduled integrally. We call the resulting problem INTEGER BOUNDED 
Ratio (IBR). In order to obtain the EPTAS for IBR, we note that in our algorithm for BR, each machine 
receives at most two small jobs fractionally. For the maximization problem of IBR we simply remove 
the fractional parts. This decreases the completion time of each machine by at most 2eWi, and thus the 
completion time of each machine is decreased by a multiplicative factor of at most 1 — 2e. This gives an 
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EPTAS for the maximization problem of IBR. For the minimization problem we assign (integrally) each 
small job to the first machine which gets a fraction of the job in the solution to BR. This may increase 
the completion time of a machine by at most eWj. Therefore, the completion time of each machine 
increases by a multiplicative factor of at most 1 + e (it may decrease for some machines as well). Hence, 
in this case the total cost of the resulting solution to IBR is at most (1 + s)p times the cost of the solution 
to BR. Thus, this gives an EPTAS for the minimization problem of IBR as well, and the following result 
is established. 

Theorem 14 Problem IBR has an EPTAS. 



4 Applying the shifting technique 

In this section we use the shifting technique of Hochbaum and Maas Il22l . We modify the total weight 
of a machine if it lies within illegal intervals. We will choose the set of illegal intervals (which imply a 
set of valid intervals) such that the following properties hold: First, the value of the objective function 
of an optimal solution using the modified weights will be close to its value according to the objective 
function of the original instance. Second, the ratio of weights of two values which are not separated by 
an illegal interval that is, belong to one valid interval will be bounded by a function 7(e). Finally, the 
ratio between the two extreme points of an illegal interval is at least ^. The set of illegal intervals is 
finite, and each such interval is bounded. 

Given such a set of illegal intervals S = {(ao, 60)) {o-i, ^i), • • • > {o-r, K)} where bi < a^+i for all 
< i < r — I and oq > min^ pj and br < Yl^=i Pj- consider a schedule and assume that when a 

machine i gets an allocation of a total weight Wi, it will contribute ^ ^^^'^ ^ to the objective function. 
We next define the function g. 

For the maximization problem where p < 1, we define g{x) = x if x ^ {ai,be) for all i, and 
otherwise g{x) = 0. For the minimization problem where p > 1, we define g{x) = x if x ^ (a^, bi) for 
all £, and otherwise if 2; G (a^, bi) we let g{x) = 2b£. By using the value of g instead of the total weight, 
the objective function value of any feasible solution may become worse (larger for the minimization 
problem and smaller for the maximization problem). We denote by p^nm the minimum size of a job in 
the instance. 

Our algorithm will choose the best outcome among a constant number of iterations. In each such 
iteration we will use a different set of illegal intervals. Let p = ^2 i l]+i ^ ^ (where equality holds for 
p < 1), then for 7? = 0, 1, . . . , p — 1, in iteration -q we will use the following set of illegal intervals: 
{a£, bi) = (pmin • (7)''^^'' )Pmin ' {j)^^^^^^), for the non-ncgativc values of i such that bi < J2^=iPj- 
The number of non-negative powers of ^ for which we define an interval in all iterations is at most 

and this is a polynomial in n, m and i. We denote by g^j the function g in the rj-th iteration, and we 
denote by Srj the set of illegal intervals of iteration t]. 

We next show that there is a value of t] such that the objective function value of the optimal solution 
to the problem with the modified weights is within a factor of 1 + e of the objective function value of an 
optimal solution to the original instance. 

Lemma 15 Denote by OPT^ the objective function value of an optimal solution with respect to the 
modified weights g^. Then, for the maximization problem there is a value ofr] such that OPT > OPT^ > 
OPT • (1 — e^) > and for the minimization problem there is a value of rj such that OPT < OPT^ < 
0Pt(1 + e). 
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We next show that given an instance where the value or cost of a solution is computed using g^, we 
can restrict ourselves to solutions which (almost) do not use weights from Sjj. More precisely we show 
the following claim. 

Lemma 16 Consider the optimal solution OPT for the instance with the modified weights g^j. Without 
loss of generality, there is at most one machine i whose weight Wi belongs to one of the intervals of Sn- 

Note that the number of machines of weight zero can be arbitrary. 

Proof. Denote by the set of machines whose weight belong to an interval from iS^. Consider 
the maximization problem. If the claim does not hold, then by assumption |M^| > 2. We modify 

the solution so that all the jobs which were assigned to one of the machines in are assigned to the 
maximum index machine of M^. The machines of did not contribute a positive value to the objective 
function in OPT^ (since the modified objective function is used), and hence this modification did not hurt 
the optimality of OPT^. The claim holds for the new solution. 

Consider the minimization problem. We repeat the following modification of OPT^ as long as 
\M^\ > 2. Let i e be a most loaded machine in OPT,, among the machines of M^. Let i' ^ i 
be another member of M^. Then, the weight of i' in OPT^ is at most Wi. We modify OPT,, by assigning 
to machine i all jobs which were assigned previously to either i or i'. This modification will at most 
double the weight of i, and hence the solution remains optimal (either i remains in and in this case 
it pays the same as in the original solution and i' pays nothing, or i is removed from and in this case 
its work is at most double the infimum point of the following allowed interval so it now pays at most its 
payment in the original solution, and i' pays nothing). Repeating the process decreases the cardinality 
of {i' is removed from after its weight becomes zero). ■ 

Next, given such a value of r], we can guess the parity of the index £ (of the interval (a^, h^) G Sr^ 
which contains a weight of a machine in OPT^ if it exists). By the previous lemma there is at most one 
such value of i. Given such a guess, we allow the use of intervals of the same parity as ^. That is, if ^ is 
an odd number, we remove from 5,, the intervals of the form (a2i-i, b2i-i), and if i is an even number, 
we remove from 5^ the intervals of the form (a2i, 624)- We denote the resulting set of illegal intervals 
by Sn^ff, where (j) G {odd,even}. The number of possibilities for this guess (of r] and (p) is polynomial 
in p Hence, we can assume that the set Sr^^^ satisfies that OPT^ does not have a machine whose weight 
belongs to an interval from 5^^^. Moreover, the following observation holds. 

Observation 17 The ratio between two weights W < W which are not separated by an interval from 

Sri^(P is bounded by a function 7(e) = ^^p_i (i.e., ^ < "/{e)), and the ratio between W' and W if they 
are separated by an interval from S^^^ is at least ^. 

Using this set of intervals 5*^, we can evaluate all the solutions which satisfy the property that they do 
not use weights from the set 5^,^, according to the weights and not the modified weights. This cannot 
harm the quality of the solution with respect to the modified objective functions. We conclude that using 
such solutions which do not use a weight from 5^ ^ does not hurt the approximation algorithm too much. 
That is, we estabUsh the following lemma. 

Lemma 18 There is a value of rj and (p, such that OPT^ does not use a weight which belong to an 

interval o/S',,^^, and such that the objective function value of OPT (with respect to the original objective 
function) in terms of the weights ( i.e., as a solution to the original instance) is within a factor ofl + eof 
OPT. 

5 Dynamic programming to approximate OPT^ 

Given fixed values of 77 and (j), the set of illegal weights S^^^ leaves a set of valid intervals whose pos- 
sible weights denoted by = {[wq, wi], [uj2, W3], . . . , w^']} where the sequence is monotone 
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increasing, ujq > Pmin, and r' < r. Recall that some machines may get a zero weight in an optimal 
solution, but we can guess their number and remove the set of these machines (of lowest indices) from 
the instance. Hence, without loss of generality it suffices to consider solutions which assign at least one 
job to each machine, and therefore the minimum weight of a machine is at least loq. 

By Observation [TTl we conclude that for any value of ?, "^^^ = - and !^^21±i < 7(e). We define a 

linear order on the intervals of U saying that an interval [uj^, tJ^+i] is smaller than [cj^/, a;^/+i] if ^ < 
and an interval [uj^,uj^^i] is at most an interval [tj^/, a;^/_|_i] if ^ < We next describe an allocation 
of jobs to intervals of Q. A job j of size pj is associated with an interval [uj^-i,uj^] if pj < uj^ and 
Pj > wg_2 where we use the convention a;_i = 0. We define an assignment of intervals of (consecutive) 
machines to intervals of O in the following sense. An interval of machines [i, i'] with parameters A, B 
is assigned to an interval [a;g,a;^+i] € if the following four conditions hold: 1) the weight of each 
machine i € [i, i' — 1] is in the interval [w^, w^+i]; 2) no other machine has weight in this interval; 3) the 
total size of jobs associated with smaller intervals and are scheduled by a machine of index at least i is 
B; and similarly 4) the total size of jobs associated with intervals at most [wg,a;^_|_i] and are scheduled 
by a machine of index at least i' is A. Here A and B are integer multiples of ^. We say that A and B are 
the parameters of the interval of machines [i, i']. 

Claim 19 The number of possibilities of a machine interval and a pair of parameters is polynomial in 
the input size. 

Note that given an interval of machines [i, i'] and values A, B we get an instance of the IBR problem 
for which we presented an EPTAS (this IBR instance has an empty set of machines if i = i' and 
otherwise at least one machine). Thus our scheme applies this EPTAS for each possibility of a machine 
interval [i, i'] such that i < i' corresponding to weight interval [w^, w^+i] with values A, B. We denote 
by IBReptas{h i' — 1, ^) B) the solution returned by the EPTAS for the IBR instance as well as its 
objective function value. If i = i' and A = B, then IBR^ptasih i' — 1, A, B) = 0, and if i = i' and 
A^ B then IBReptasih i' — B) returns FALSE. If the returned output is FALSE, then the value 

is 00 for the minimization problem and —00 for the maximization problem. 

To find the approximated solution for the minimization problem we find a shortest (minimum cost) 
path in the following layered graph G, and to find the approximated solution for the maximization 
problem we find a longest (maximum cost) path in this graph G (since G is a layered graph, it is acyclic 
and hence both the shortest path problem and the longest path problem are solvable in linear time). 
We have a layer for each value of ^ such that [a;g,a;^_|_i] E Q. (that is, we will have a layer for every 
even value of 0. Each such layer corresponding to ^ has a vertex for each machine (and one additional 
dummy machine of index m + 1), and each possibility for the value of A (i.e., for each integer multiple 
of Given a vertex (i, b) in layer ^ and a vertex {i' , a) in layer ^ + 2, there is an arc from the former 
vertex to the later vertex ifi< i', and the cost associated with such an arc is IBR(.ptas{i, i' — 1, 6). 
The construction of G takes polynomial time, and it has a polynomial size. 

We next find a shortest or longest path in G from the vertex (1,0) of the layer with index to the 
vertex (m + 1, 0) in the last layer. We schedule the jobs according to the path which we found. That 
is, if the path uses the arc from {i, b) of layer ^ to {i', a) of layer + 2 where i < i', then we schedule 
large and small jobs as defined by the solution IBR^ptasih i' — ^, a, b). The total size of jobs which 
are associated with the interval [tj^, (xj^+i] or smaller intervals, and are not scheduled to machines of 
index at most i' — 1 is indeed at most a for the minimization problem and at least a for the maximization 
problem. The EPTAS for IBR may use at most one small job which exceeds the total of B. This job is 
scheduled in the case of minimization to the machine i' — 1 and discarded in the case of maximization. 
In addition, since we use only some values for w{K), in the case of minimization we have too large 
room for small jobs, and in the case of maximization too small room for small jobs. Thus, for the 
minimization problem, it can be the case that too many jobs are scheduled and as a result no small 
available jobs remain. In this case, the algorithm moves to deal with the next arc of the path (or stops 
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if this is the last arc). Similarly, for the maximization problem, since discarded jobs remain available, 
there are always sufficiently many small jobs to cover all the configurations (and some small jobs may 
remain unassign upon the termination (at the end of the path). Note that every job is scheduled by one 
of the solutions corresponding to this path, and hence we obtained a feasible solution whose objective 
function value is the total cost of the arcs in the path. 

We next note that OPT^ also corresponds to a path in G in the following sense. If OPT^ uses machines 
i, i + 1, . . . , — 1 with weight in the interval [w^, w^+i], the total size of the jobs allocated to smaller 
intervals and are not scheduled by OPT^ to machines with index at most i — 1 is 6, and the total size 
of the jobs allocated to intervals at most this interval and are not scheduled by OPT^ to machines with 
index at most i' — 1 is a, then we say that the arc from (i, b) of layer ^ to (?;', a) of layer ^ + 2 belongs 
to the path associated with OPT,,. Since OPT^ is a feasible solution, this set of arcs can be augmented to 
form a path by adding zero cost arcs (from (i, b) of layer ^ to (z, b) of layer ^ + 2). Moreover, the cost 
of OPT^ can be also partitioned into its arcs, by assigning the cost (according to OPT^) of the machines 
i. i + 1, — 1 to the arc from (i, b) of layer ^ to a) of layer ^ + 2. By the correctness of the 
EPTAS to IBR we conclude that the cost assigned to such an arc is within a multiplicative factor of 1 + e 
of the cost of the arc. Hence, the cost of this path in G is within a factor of 1 + e of the objective function 
value of OPT^. Since we use the optimal path, we are not worse than this path of OPT^, and thus we 
established the following result. 

Theorem 20 Both the problem of minimizing ^ Cf and the problem of maximizing Yl^i for 
real finite values ofp have efficient polynomial time approximation schemes. 
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A Omitted proofs 
A. 1 Proof of Lemma |2] 

Assume by contradiction that Wj > 5Wii. The solution pays at least for this pair of machines. 

We move all the jobs which were scheduled on machine i to machine i' . It suffices to show that the 
completion time of machine i' in this new solution is smaller than This last claim holds because the 

completion time of i' in the new solution is — < < ■a{5) = where the first inequality 

holds because Wii < ^ and 5 < 1, the second inequality holds because Sj < a{6)- Sj', and the equality 
holds by the definition of a{6). 



A.2 Proof of Lemma |3] 

Assume by contradiction that the claim does not hold, that is, there exist such machines i and i' with 
Wi > 5Wii . We compare the current solution with a new solution that schedules all the jobs (which 
were previously scheduled on either i or i') to machine i' . It suffices to show that this new solution 
is better, that is, that the following inequality holds: ( j > | j _|_ [ ) . We first note 



that < (5t)^ • iC^ + ^y - !)■ This claim holds because Wi < Wi' and Si> < a{6) ■ Si, 

. . . /tj/ \P /W , \P , _ . . . _ 



and therefore ^ ('^('^))^' ^^'^ claim holds by the definition of a{d). Therefore, 

to get a contradiction it suffices to show that f ^'-^^i' \ > f ' (1 + which is equivalent 



to {Wi + Wi'Y > (Wi' ■ (1 + 6)y which holds using Wi > 5Wi', since is a strictly monotone 
increasing function of x (for < x < oo and any p > 0). 



A.3 Proof of Lemma 13 

The first claim holds because we round down the size of each job. We prove the second claim. Consider 
a machine i whose work in SOL is denoted by Wi. If Wi < np then its new work is at least 0. Otherwise, 
its new work is at least Wi — nfi. Consider the case where Wi > nfi. Then, the contribution of machine 

^ to SOL_ is at least [m^Y > {0 - (f = (f )' " > (f )' " where the 

first inequality holds by concavity of x^' for p < 1, and the last one holds using Si > Sm- 

Consider the case where Wi < nfi. Then, the contribution of machine i to SOl^new is at least 
> ~ (i^) ^ (^) ~ ' where the second inequality holds similarly to the previous 

case. Thus, the decrease of the contribution of each machine towards the objective function value of 
SOL, is at most < where the inequahty holds because SOLow > {^f^^ which is the 

value of a solution that places the largest job on the slowest machine, and ignores all the other jobs. 



A.4 Proof of Claim [To] 

First, note that OPT^r induces a feasible solution for the integer program. To get this induced solution 
we apply the following rounding for each machine in the optimal solution: for each machine we define 
a configuration K by first counting the number of large jobs of each size, and it remains to define w{K). 
w{K) is a rounded up value of the total weight of the jobs assigned to the machine to the closest integer 
power of 1 + e. We next argue that the total work of this machine in OPT^r is at most w{K). This is 
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so because otherwise, a machine with total work of 6 in OPT{,r has a configuration K with w{K) > 6 
which is an integer power of 1 + e, and w{K) = 9' such that 6' < 9. However, since 9 is an integer 
multiple of fi we have 6* = /i • [^J < ^ - [^^\ = 9', that is a contradiction. The total size of the gaps 
for scheduling the small jobs does not decrease, and thus the total size of the unscheduled jobs is at most 
A. This results in a feasible solution to the integer program. The value of the objective function of this 
rounded solution is at most (1 + e)P ■ OPTft^, and at least X*, and thus the claim holds. 

A.5 Proof of Claim [U 

First, note that OPT^r induces a feasible solution for the integer program. To get this induced solution 
we apply the following rounding for each machine in the optimal solution: for each machine we define 
a configuration K by first counting the number of large jobs of each size, and it remains to define w{K). 
w{K) is a rounded down value of the total weight of the jobs assigned to the machine to the closest 
integer power of 1 + e. We next argue that the total work of this machine in OPT;,r is at least w{K). This 
is so because otherwise, a machine with total work of 9 in OPT;,,, has a configuration K with w{K) < 9 
which is an integer power of 1 + e, and w{K) = 9' such that 9' > 9. However, since 9 is an integer 
multiple of we have 6* = /i • [^] > fi ■ l^^] = 9', that is a contradiction. The total size of the gaps 
for scheduling the small jobs does not increase, and thus there is sufficient total size of small jobs to fill 
in all the gaps, and the total size of the unscheduled jobs is at least A. This results a feasible solution to 

OPT 

the integer program. The value of the objective function of this rounded solution is at least jjif^, and 
at most X* , and thus the claim holds. 

A.6 Proof of Claim [12] 

First, the construction of the integer programs takes polynomial time since the set of all configurations 
can be enumerated in polynomial time (using the fact that JC has at most a constant number of configura- 
tions which is upper bounded by a function of e). Next, we observe that the dimension of each of these 
programs (the number of variables) is |/C| + \H\, and as explained above, both |/C| and | are bounded 
by a function of i. Thus the integer program has a fixed dimension, and we can use the polynomial time 
algorithms for solving such a problem. The number of constraints (beside the non-negativity constraints) 
is r(e) + \H\ + 2 which is again bounded by a function of ^. Therefore, using Lenstra's algorithm \2T} 
or one of its improvements, give a polynomial time algorithm for solving (exactly) the integer programs 
(recall that the time complexity of solving an integer program of dimension d is f{d) - poly where / is an 
exponential function of the dimension, and poly is a polynomial in the binary encoding of the program). 
To obtain a strongly polynomial time we use the following observations. First, the coefficients in the 
objective function are integer powers of 1 + e and can be scaled to be at most a function of e (since the 
ratios and ^^t/~^ are bounded). Next, we scale constraints (|3]l, (fUl, and (|9ll by dividing the 
constraints by the factor /i. In the resulting constraint matrices and right hand sides, all the coefficients 
are strongly polynomial (i.e., do not depend on the magnitude of the numbers in the instance). 

A.7 Proof of Lemma [15] 

Fix an optimal solution OPT to the original instance, and associate the value of each machine to a 
corresponding value of rj as follows. Assume that a machine i has weight Wi, and assume that Wi € 
(Pmin ■ (i)^"'^^'' ,Pmm ' (1^^"^^''^^), then wc associatc (^^^ with this value 7]q of r]. Note that for a 
machine i there is at most one such corresponding value of rj. If there is no such corresponding value of 
r] then we associate (^^^ with an arbitrary value of rj. By the pigeonhole principle, there is a value r/o 
of r? which is associated with at most 



14 



Consider first tlie maximization problem. Using the value rjo of rj, OPT is a feasible solution to the 
problem with the modified value function gr/g whose objective function value is at least OPT • (1 — = 
0Pt(1 — e'^), and thus the claim holds. 

Next, consider the minimization problem. Using the value 7]q of t], OPT is a feasible solution to 
the problem with the modified cost g^fg whose objective function value is at most OPT + ■ (|)^ < 
Qprp^ 0£T . _i_ ^ 0PT(l + e), where the first inequality holds because the modified weight of a machine 
whose corresponding value of r] is r/o is at most | times its weight. Thus, the claim holds. 

A.8 Proof of Claim [19] 

The number of pairs of machines i and i' is at most m?, and the number of possibilities for the value A 
(and similarly for B) is at most '-^^^^^ < ■ Thus, the claim holds. 
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